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Abstract 

We perform a systematic study of the effects of the type-I seesaw mechanism on the dark 
matter abundance in the constrained super symmetric standard model (CMSSM) which in- 
cludes three right-handed neutrinos (the i^CMSSM). For large values of mQ,mi/2, we exploit 
the effects of large neutrino Yukawa couplings on the renormalization group (RG) evolution 
of the up-type Higgs. In particular, we show that the focus point scale can greatly exceed the 
electroweak scale resulting in the absence of a focus point region for which the relic density 
of neutralinos is within the range determined by WMAP. We also discuss the effects of the 
right-handed neutrinos on the so-called funnel region, where the relic density is controlled by 
s-channel annihilations through a heavy Higgs. For small values of mo,w.x/2) we discuss the 
possibility of sneutrino coannihilation regions with an emphasis on the suppression of the left- 
handed slepton doublet masses due to the neutrino Yukawa coupling. We consider two types 
of toy models consistent with either the normal or inverted hierarchy of neutrino masses. 

PACS: 12.60. Jv 



1 Introduction 

The seesaw mechanism offers a compelling explanation for the tiny masses of the left-handed neutri- 
nos by introducing heavy gauge-singlet right-handed neutrinos which can alleviate the fine-tunings of 
the neutrino Yukawa couplings [1]. If the right-handed neutrino masses are close to the grand unifica- 
tion (GUT) scale, one typically finds neutrino Yukawa couplings of order unity. These heavy degrees 
of freedom may also influence such effects as lepton flavor violation [2-10] and the mass/coupling 
spectrum due to changes in the RG (renormalization group) evolution [11-18]. In this paper, we 
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implement the seesaw mechanism in the context of the constrained minimal supersymmetric model 
(CMSSM) by adding three right-handed neutrinos Ni. We study the consequences of the near- 
GUT-scale right-handed neutrinos in such a model (hereafter called the i^CMSSM) on the thermal 
dark matter relic abundance which is calculated around the electroweak scale. The effects of the 
seesaw mechanisms on the dark matter abundance has been studied for several GUT and seesaw 
schemes [19-24] indicating signatures different from those in conventional CMSSM scenarios [25-32], 
such as the sneutrino coannihilation regions [18]. 

Here, we present a systematic study of the effects of a general type-I seesaw without any ad- 
ditional constraints from a specific GUT model which may require Yukawa coupling unification. 
We show in detail the behavior of two of the standard CMSSM regions where the WMAP relic 
density [33] is achieved, namely the funnel region, where the mass of the lightest supersymmetric 
particle (LSP), the neutralino, is roughly half the mass of the pseudo-scalar Higgs boson [25,27], 
and the focus point region which borders the region where radiative electroweak symmetry breaking 
is no longer possible [28,29,32]. In the case of the latter, the behavior of the Higgs mass drastically 
changes in the z/CMSSM because of the change in the focus point scale from the RG evolution of 
m^^. As we will see, the focus point region can even completely disappear once the energy scale of 
the focus point is sufficiently greater than the weak scale. As a result the WMAP strip associated 
with the focus point can cease to exist once the right-handed neutrino mass scale, M^r, is sufficiently 
large. The funnel region is also affected (though not as dramatically) as the pseudo-scalar mass is 
also affected by the choice of M^. 

We also detail changes in the stau coannihilation region which occurs in the CMSSM near the 
line of degeneracy between the neutralino and lighter stau masses. As we have shown previously [18], 
the introduction of a right handed neutrino can affect the running of the left-handed sleptons so 
that they become lighter than the more common next-to-LSP, the predominantly right-handed stau. 
This has the effect that the sneutrino may well become the NLSP. Here, we also extend that work 
to include the effects of all three right-handed neutrino states. 

We review the essential features of the i/CMSSM in ^ In ^ we consider the effect of the 
right-handed neutrinos on the location of the focus point regions in the i^CMSSM. We preform an 
analogous analysis for the Higgs funnel regions in ^ The sneutrino coannihilation regions with 
small mo,mi/2 are studied in ^ followed by the conclusions in §H1 

2 z/CMSSM 

For simplicity, we do not consider flavor mixings in the neutrino sector, and as such, in the basis 
where the heavy right-handed neutrino mass matrix is diagona0 our z/CMSSM superpotential can 
be written as 

W = WcMSSM + yNML^Hu + \mn^N^N^. (1) 

^We can always make the right-handed neutrino mass and lepton Yukawa couphng matrices diagonal by unitary 
transformations of TV^, L and E'^ . The diagonal form of the neutrino Yukawa coupling matrix in such a basis is the 
simplification we make here. 
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We assume that all soft SUSY breaking parameters including the soft SUSY breaking sneutrino 
masses, m^Vi, are universal at the GUT scale with value mo- In addition, we also assume the SUSY 
breaking trilinear couplings including A^. are universal at the GUT scale. As in the CMSSM, 
we take universal gaugino masses. Then, in addition to the conventional CMSSM parameter set 
mo, 1711/2, Ao, tan P, sgn{^), the i/CMSSM has the following additional parameters 

Mni (QgUt) , MjV2 (QgUt): Mn^ (QgUt) , m^i {Qmz ) > '^1.2 {Qmz ) : {Qmz ): (2) 

where the heavy right-handed neutrino masses M^. are specified at the GUT scale Qgut, while the 
light left-handed neutrino masses are specified at the weak scale, Qmz- These boundary conditions 
are sufficient to determine the neutrino Yukawa couplings y^^ as well. In this analysis, we restrict 
the right-handed neutrino masses to lie below the GUT scale Mjv^ < Qgut- We evolve the full two- 
loop renormalization group equations (RGEs) in the z/CMSSM. Each of the right-handed neutrinos 
decouples at Mjv^ which itself runs from the GUT scale to M^. . Below the energy scale Mj^. , is 
integrated out and the effective Lagrangian includes the residual higher order operators suppressed 
by M^Vj. Because we are interested in the value of Mj^- relatively close to the GUT scale in the 
type-I seesaw scheme, those non-renormalizable operators do not affect the mass spectrum at the 
low energy scale other than the left-handed neutrino masses which receive a dominant contribution 
from the dimension five operator 

Ls 9 -k{LHu){LH^) (3) 

Hence, we keep the running of k in our full two-loop RGEs and K{Hy)'^ is matched to m^{Qmz) by 
our electroweak scale boundary conditions [34, 35] . 

Even though our numerical analysis was performed including full two-loop RGEs, it will be 
useful to use the one-loop RGE expressions for a qualitative discussion. The relevant RGEs are 
{t = logQ) 
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(13 
Tr SyuyI + yjvyJv + ^yIYu + YdYd - —gf - 3^2 - y^s J + - 

S = Tr{xnl + nil-2nil-nil + inl)+ml^-ml^ (4) 



where, for the models to be discussed below, hij = AijUijdij. In the CMSSM as well as in the 
z/CMSSM, S" = at the GUT scale due to the universality of the soft scalar masses. Deviations 
from 5 = are due to the RG evolution at the two-loop level, hence, S does not play a significant 
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role in our study. In the following sections, the RG effects of large neutrino Yukawa couplings on 
mjj^ are studied in the focus point and Higgs funnel regions whereas those effects on the slepton 
doublets m| are studied in the sneutrino/stau coannihilation regions. 

We leave the question of flavor mixings and the CP violating phases for future work and we 
present here the analysis for two simple toy models to illustrate the effects of the heavy right- 
handed neutrinos on the thermal dark matter relic abundance. The first model includes only 
one right-handed neutrino in the third generation A^3 consistent with a normal mass hierarchy 
m^^ S> m^^ = m^^. To see the effects of the additional generations with multiple right-handed 
neutrinos, our second model includes two heavy right-handed neutrinos in the first two generations 
A^i , A'2 leading to an inverted mass hierarchy spectrum m,^.^ -C m^-^ = niy^ . 

We start our discussion with the first model including only and and concentrate on large 
values of mo,mi/2 in sections ^ and ^ In section ^ we go on to discuss our second model 
including two right-handed neutrinos in the first two generations A^^i, A^2 at relatively small mo, mii2 
relevant for the sneutrino/stau coannihilation regions. 



3 The focus point region 

When the Higgsino mass parameter, /x, becomes small or comparable to mi/2, the lightest neutralino 
(which is typically the LSP) has an increasingly large Higgsino component and neutralino annihi- 
lations become dominated by W^W~ final states through t-channel chargino exchange. The result 
of the enhanced annihilation cross section is a thermal neutralino relic abundance in the desired 
WMAP (95% CL) range [33] 

0.0975 < VL^K^ < 0.1223 (5) 

in a region (large mo) where nominally the relic density is expected to be large. This region is close 
to the limit (at large mo) where there cease to be solutions to the Higgs potential minimization 
conditions corresponding to radiative electroweak symmetry breaking. In the CMSSM, the RGEs 
exhibit a focusing property leading to what is often referred to as the focus point or hyperbolic 
branch regions [28,29,32]. The RG trajectories of rnj^^ ~ — for different values of mo (and 
tan/3) meet (or focus) near the weak scale [28]. The focusing of the RG trajectories of m'jj^ also 
occurs in the z/CMSSM, as the RGEs for the gaugino masses and trilinear A parameters remain 
independent of the scalar masses at the one-loop level (while those for the scalar masses depend on 
the scalar masses, gaugino masses and A parameters). In the i/CMSSM, however, the focus point 
scale where RG trajectories meet is not necessarily found at the weak scale due to large neutrino 
Yukawa contributions. The RG evolutions of rriHu the CMSSM and i/CMSSM are illustrated in 
Fig. [T]for different values of mo, where the sign of rriHu indicates the sign of rnfj^. We use a top mass 
rrit = 173.1 GeV [36] and a bottom mass in the minimal subtraction (MS) scheme mhlrrif,)^^^ = 4.25 
GeV in our analysisl^ 



^See, for instance, Refs. [27,28,37,38] for the effects of the variations of nit, mi, within the experimentally allowed 
ranges on the dark matter abundances, which however do not affect our discussions in this paper. 
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Figure 1: The RG evolutions of with different values of ttiq for the CMSSM (dashed) and 
the z/CMSSM with Mj^^ = lO^^GeV and m^^ = O.OSeV (sohd). We have set mi/2 = 300 GeV, 
tan P = 10, Aq = 0, and /i > 0. 

Before presenting the numerical results for the dark matter abundance calculations, we briefly 
examine the RG evolution of mjj^ to illustrate the properties of the focus point regions which yield 
the observed thermal relic abundance. The tree-level minimization of the Higgs potential gives 

" tan2/3-l 2 ^ ' 

This allows one to estimate /i^ ~ ""^l/u ~ moderate to large values of tan (3 (we hereafter 

assume fi > in our analysis as favored by measurements of 6 — 57 and the anomalous magnetic 
moment of the muon). The RGE of mfj^, which gives us an estimate for fi at the weak scale, is given 
in Eq. (jlj). In the CMSSM (or when = 0), we see the RG evolution of m|^^ is sensitive to the top 
Yukawa coupling [39,40]. In fact, in the CMSSM, the focus point scale Qfp is determined by the 
value of the top quark Yukawa coupling evaluated at the focus point and the (entire RG trajectories 
of) gauge couplings [28]. The focus point scale of m^^ in the CMSSM turns out to be the weak 
scale for the measured value of top quark mass with a small dependence on moderate to large values 
of tan/3 (this can be inferred from the top Yukawa coupling yt oc 1/ sin/3 oc const + 0(tan~^/3)). 
This behavior is seen in Fig. [1] which shows the CMSSM running of mfj for three values of mo. 
As one can clearly see, the three dashed curves corresponding to the CMSSM focus at a weak scale 
value for Q. 

In the t/CMSSM, the effects of the neutrino Yukawa couplings on the RGEs of and yt can 
be significant as seen in Eq. (jl]), so that the focus point scale of m'jj^ can well exceed the weak 
scale. This is seen in Fig. [1] which shows the three solid curves of the z/CMSSM focusing at the 
scale ~ 10^ GeV where yNaiQcur) ~ 2. Consequently, and contrary to CMSSM, the insensitivity of 
m|^^ evaluated at the weak scale to uiq and tan P does not necessarily hold in the z/CMSSmH The 

■^Once the heavy right-handed neutrino is integrated out at Q = Mn, the analytical estimation of the focus point 
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positions and existence of the focus point regions (at large mo where the rehc density is sufficiently 
small) change due to this shift in the focus point scale. When the focus point scale is below the 
energy scale relevant for the neutralino annihilations at which n is evaluated, we can decrease 



II ~ \J —^"h^ by choosing a larger value for mo (see CMSSM curves in Fig. [T]). Recall that at small 
/i, neutralino annihilation is dominated by the Higgsino component and can yield an acceptable 
relic density. As one can also see from Fig. [H when the focus point scale becomes large (as depicted 
in the figure for the i^CMSSM), the weak scale value of Tn'\j^ becomes sensitive to mo and is driven 
to large negative values when mo is large. In this case an increase of mo (with other parameters 
fixed) increases and the Higgsino component of the LSP is diminished and no region of suitable 
relic density can be found. 

We also note that, for a given value of mo, the value of mfj^iQ) is smaller in the z/CMSSM 
than in the CMSSM, as can be inferred from the form of RGE of m'jj^ in Eq. (jlj) due to the 
additional terms involving yj^. The behavior of these mass parameters is also illustrated in Fig. [2] 
where we also show the shift/disappearance of the boundary where radiative electroweak symmetry 
is no longer possible. Inside the (pink) shaded region, the Higgs potential minimization conditions 
yield /i^ < 0. Along the boundary of that region /i = 0, and close to the boundary, we show a 
contour of constant n = 200 GeV. Almost coinciding with that contour is the very thin WMAP 
focus point region (shaded turquoise) where the relic density lies in the range given in Eq. ([5]). Also 
shown are contours of constant /i at 200 GeV intervals. As one can see, as M^^ is increased, the 
focus point region shifts to increasing mo and is absent at M^.^ ~ 10^^ GeV (which corresponds to 
UNaiQcur) ~ 0.6) at which point there is no longer any sensitivity of the value of /i to mo. This 
corresponds to the coincidence of the focus point scale with the scale at which /i is evaluated. In 
this figure, we have fixed mi/2 = 300 GeV, tan/3 = 10 and Aq = 0. For M^^ ^ 10^^ GeV, the focus 
point energy scale, Qfp is low and increasing mo helps making m|^^ less negative which reduces /z 
allowing for the presence of the focus point region. For larger M^vg, the opposite behavior occurs 
and increasing mo increases /z. 

Fig. [2] also shows the region where the two-loop contributions become large. To the right of 
the vertical dashed line at M^.^ ~ 2 x 10^^ GeV, yl^^/lGvr^ > 1/10 (where the dominant two loop 
contribution is roughly 1/4 of the 1-loop contribution). At larger values of M^s, our perturbative 
calculation results could become unreliable [4,7,41]. To the left of this line, we see that /i increases 
with M7V3 for fixed mo. Finally, in the region below the (red) dot-dashed curve, rrih < 114.4 GeV and 
fails the LEP limit [42] (modulo the uncertainty in the calculation of rrih for which FeynHiggs [43] 
was used here). 

Our qualitative discussion so far can be summarized in Fig. [3] which shows the (mo,mi/2) plane 
for tan/5 = 10 and Aq = 0. In the lower right corner of the figure, the brown shaded region 
corresponds to parameters for which the LSP is the charged partner of the tau lepton and as such 
is excluded. Running along this region is the so-called WMAP co-annihilation strip, where co- 
annihilations between neutralinos and staus are largely responsible for obtaining the WMAP relic 

analogous to CMSSM can be applied to the z^CMSSM for Q < with non-universal scalar mass values at Q = M^. 
Hence, as in the case of the CMSSM with non-universal boundary conditions, the focus point can still be found at 
the weak scale in the z/CMSSM if one can tune the model parameters to realize the particular ratios of those scalar 
mass values at Q = Mn [28]. 
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Figure 2: The (mojMTVg) plane for fixed mi/2 = 300 GeV, tan/3 = 10, = 0, and m^g = 0.05 eV. 
Tlie (pink) shaded region shows the portion of parameter space where the radiative electroweak 
symmetry breaking conditions can not be satisfied. The focus point region (turquoise) compatible 
with Eq. ([5]) runs along this boundary nearly coinciding with the contour for fi = 200 GeV. Other 
contours of constant /i are also shown. The (red) dot-dashed curve shows the contour of rrih = 114.4 
GeV. Below this curve, the Higgs mass is below the LEP limit. To the right of the vertical dashed 
line, y^g/lGvr^ > 1/10 and our perturbative expansion becomes suspect. 

density. This area (as in the previous figure) is shaded turquoise. In the lower left corner of the 
figure, the small dark (green) shaded region is excluded as the supersymmetric contributions to 
6 — s> S7 disagree with the experimental determination [44], while the light pink shaded region is 
favored by the measurement of the muon anomalous magnetic moment at the 2-a level [45]j. The 
region between the pair of dashed line gives agreement to within la of the experimental result. The 
red dot-dashed contour corresponds to a Higgs mass of 114.4 GeV. At lower mi/2, the Higgs boson 
would be lighter, which is excluded by its non-observation at LEP [42] . We also plot a black dashed 
contour for m^± = 104 GeV, the region at lower mi/2 also being excluded by LEP. These contours 
are shown for four choices of Mjvg as labelled. Note that the stau co-annihilation strip is barely 
sensitive to the choices of M7V3 made here as are other quantities at low rriQ. This is reasonable 
because the terms involving a neutrino Yukawa coupling in the RGE of the slepton doublets can 
become large when the scalar masses or trilinear A terms are large as seen in Eq. (jlj). Because the 
gaugino masses and the right-handed stau are not affected by a neutrino Yukawa coupling at the 
one loop-level, the ratio of the stau and gaugino mass would not be affected significantly by the 

''We have used (5a^ = 24.6 ± 8.0 for the discrepancy between theory and experiment based on e+e~ data [46]. 
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right-handed neutrino when mo and A are small. 

In the upper left corner of Fig. [3l we see the (pink) shaded region where there is no radiative 
electroweak symmetry breaking. Along the boundary of this region, we find the focus point strip 
where again the relic density of neutralinos is within the WMAP range. The shaded region corre- 
sponds to the standard CMSSM model or equivalently the i^CMSSM with Mn-, = Qgvt ^ 2 x 10^*^ 
GeV. Within the shaded region, we also show how the region of no radiative electroweak symmetry 
breaking recedes to higher rriQ for values of M^rg < Qgut- The boundary in each case is shown 
labelled by the chosen value of Mjvg . The WMAP focus point strip would run along each boundary. 
Note that for very low M^.^ (~ 2 x 10^^ GeV), we effectively recover the CMSSM as the neutrino 
Yukawa couplings are too small to have any effect on RGEs. As M^^ is increased, we see the 
disappearance of the focus point region as discussed above. 




Figure 3: The (mo, mi/2) plane for tan/5 = 10, Aq = 0, m^^ = O.OSeV, and four choices of M^^. The 
region labelled Qgut corresponds to the CMSSM. The black solid lines represent the boundaries of 
the regions above which (i.e. for bigger mo) there is no electroweak symmetry breaking and < 0. 
Contours and shaded regions are described in the text. 



4 Higgs funnel region 

Another efficient neutralino annihilation process which results in the requisite thermal relic density 
is through the s-channel Higgs (A and H) resonances and is possible at large tan /3. The dominant 
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contributions are through the pseudo-scalar A which is CP odd and can couple to initial s-wave 
states (while the couplings through the CP-even H suffer from p-wave suppressions due to the small 
velocity of the neutralinos) , with the dominant final states bb whose coupling is enhanced by a 
factor oc mfetan/3. At large tan/5, the resonances become broad enough (typically 10 ~ 50GeV) so 
that sufficient neutralino annihilations can occur even if is several partial widths away from the 
exact resonance region 2m^ = tua- The change in the RG evolution of m|^^ in the i/CMSSM also 
affects the A-pole funnel regions. This can be understood from the tree-level Higgs potential which 
gives 

= + "^k + 2|/xp (7) 
which, for moderate to large tan /5, becomes 

For fixed mo, and noting that the gaugino mass RGBs include i/at at the two-loop level, we would 
expect the Higgs funnel regions to move towards larger mi/2 for larger y^. This is because 
increases vha by making m|^^ more negative and hence requires larger mi/2 in order to satisfy 
~ vha- For mi/2 fixed, the increase in from an increase in t/at must be compensated by a 
decease in mo. This is because a decrease in mo can change m|^^ (which is affected by yj<! at the 
two-loop level) to cancel the change in m|^^ and can also suppress the terms in the RGE which 
involve the products of y\q and scalar mass squared. These effects are summarized in Figs. H] and 
O Fig. IHshows the contours of m^/2m^ in the (mg, M^.^ plane. Here, we can see the broad Higgs 
funnel region extending between the mA/'2,m-^ = 1 and 1.1 contours. The Higgs resonance is so 
efficient such that the WMAP strips typically show up slightly off the exact resonance myi/2m^ = 1 
contour line. In the brown shaded region at lower mg, the lightest stau is the LSP and hence this 
region is excluded. Between this region and the WMAP strip, the relic density is lower than the 
WMAP range. 

Fig. [5] shows the (mo,mi/2) plane at tan/? = 55 for which the funnel region is clearly visible. 
As in Fig. [3l we show the position of the Higgs and chargino mass contours at the LEP limit. 
Also shown are the shaded regions excluded by a charged stau LSP and 6 — S7 as well as the 
shaded region preferred by — 2. The shaded region in the upper left again corresponds to the 
region with no radiative electroweak symmetry breaking in the CMSSM. The shaded region labelled 
CMSSM, corresponds to the funnel region where the relic density agrees with the WMAP range. 
The remaining funnel-like shaded regions correspond to the shift in the funnel in the i/CMSSM for 
different values of Mn^ = 2 x 10^^ GeV, 5 x 10^^ GeV and 10^"^ GeV. Fig. H corresponded to the 
choice mi/2 = 1100 GeV, and as one can see for values of mo above the line of neutralino-stau 
degeneracy and below the WMAP funnel, s-channel neutralino annihilations via the heavy Higgs 
are so strong near ~ m^ that the relic density is pushed to values far below the WMAP range. 
As such, for this value of mi/2 there is no WMAP stau coannihilation region. 

As the neutrino Yukawa coupling becomes larger, the funnel regions move down towards the 
stau coannihilation regions. Also note that the larger values of /i obtained through the effects of 
the neutrino Yukawa couplings increase the left-right mixings in the stau mass matrix, especially 
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Figure 4: The (mo, M^^) plane for fixed mi/2 = 1100 GeV, tan/3 = 55, Aq = 0, and m,^^ = 0.05 eV. 
The funnel region (turquoise) compatible with Eq. ([5]) runs almost horizontally around mo — 1100 
GeV. Contours of constant m^/2m^ are also shown. To the right of the vertical dashed line, 
l/^g/lGvr^ > 1/10 and our perturbative expansion becomes suspect. 

for large tan/3, due to the off-diagonal term — mT-(/i tan /3 + A^-). The neutrino Yukawa couplings 
can also increase the left-handed component in the lighter slepton mass eigenstate by suppressing 
the left-handed slepton mass as discussed in §5.21 

The focus point strip is also visible in Fig. [51 For the case of the CMSSM, in the area between 
the focus point strip and the shaded region corresponding to //^ < 0, the lightest neutralino acquires 
a significant Higgsino component and the relic density is below the WMAP range. As discussed in 
the previous section, as M7V3 is increased, the focus point strips move to larger values of mo. The 
black solid lines to the left of each colored strip correspond to the edge of the region where radiative 
electroweak symmetry breaking is lost. 

5 Effects from the first two generation right-handed neu- 
trinos 

We have so far discussed the effects of only one heavy right-handed neutrino in the third generation 
exemplifying a normal hierarchy scenario m^^ ^ mj^i ^ m^2- That is we have implicitly assumed 
that the Yukawa couplings for the first two generation neutrinos are small and do not affect the spar- 
ticle spectrum. We next consider the effects of multiple right-handed neutrinos in other generations 
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Figure 5: As in Fig. [3], the funnel regions for the CMSSM and for different values of Mjvg = 
2 X 10^^ 5 X 10^3^ 10^^ GeV in the z/CMSSM. 



for comparison. We then go on to explore how the parameter space at relatively low mi/2 and mo is 
affected by the seesaw mechanism, with an emphasis on the emergence of sneutrino coannihilation 
regions. 

5.1 Focus and Higgs funnel regions 

We now consider a model which includes A^^i and N2 with equal masses, = Mjvj , corresponding 
to an inverted mass hierarchy irty^ <^ niyi = m,^2- Even though this toy model is not complete 
in that we ignore the neutrino mixings for simplicity, it does illustrate the features of the e and 
H right-handed neutrinos which are distinct from those of the r right-handed neutrino and those 
of the CMSSM. We save a more complete analysis including all three neutrinos and mixings for a 
future publication. 

The effects of Ni and N2 on the Higgs masses are qualitatively analogous to those of because 
the neutrino Yukawa coupling term, yijNiLjHu, {yij is diagonal in this example) which gives the 
key interactions between and the Higgs bosons has a common form for each generation. Thus we 
expect similar effects on the RG evolution of m|^^ (or /i) as discussed above. The effect of making 
rnjj^ even more negative is, however, bigger now because there are contributions from both first and 
second generations, and hence the shift and disappearance of the focus point and funnel regions 
are more prominent quantitatively. As one can see from Fig. [6], the edge of the region with no 
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electroweak symmetry breaking at large mo is shifted down in M^^ by a factor of two relative to 
where it was in Fig. [2] for N'^. The sensitivity of /i on M^r^, M7V2 is also shown in Fig. [61 




Figure 6: As in Fig. [2], the (mo,MArJ plane for fixed mi/2 = 300 GeV, tan/3 = 10, = 0, and 
rriu^ = rriy^ = 0.05 eV. It is assumed that M^-^ = M^^. 

5.2 Sneutrino coannihilation regions 

Our final consideration is the sneutrino coannihilation region at relatively small mo and mi/2- Here, 
the qualitative interpretation of the parameter space compatible with the observed dark matter 
abundance can differ significantly between the two toy models discussed in this paper. Sneutrino 
coannihilation regions where the sneutrino is the NLSP with a neutralino LSP are not realized 
in the CMSSM because the right-handed species are usually lighter than the corresponding left- 
handed slepton doublets. This is because the SU(2) gauge coupling terms in the RGBs tend to push 
the masses of slepton doublets up over their right-handed counter parts when universal boundary 
conditions at the high energy scale are chosen as in the CMSSM. In the z/CMSSM, however, the 
neutrino Yukawa coupling which involves the left-handed SU(2) doublet unNLHu can push the 
left-handed slepton doublet mass down while running from Q = Qgut down to Q = M^, enabling 
the slepton doublets to be lighter than their corresponding right-handed singlets at the electroweak 
scale. 

In the model with a single right handed neutrino {N^), the realization of the sneutrino coanni- 
hilation regions tends to require a large value of Aq and a moderate value of tan/3 [18]. A small 
universal scalar mass is generally required to make the sneutrino light, which results in the necessity 
for large Aq because the effects of the neutrino Yukawa coupling show up in the RGBs in the form 
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of a product between un and the sum of scalar masses and Aq as given in Eq. (jl]). Large Aq also 
helps enhance the stop loop contributions to the Higgs mass which relaxes the tight constraint from 
Higgs mass lower bound when mi/2 is small (low mi/2 is also preferred in order to obtain a light 
sneutrino by suppressing the the RG evolution of the slepton doublet due to gauge interactions). 
Moderate tan/? is preferred because at large tan/3 the induced mixing in the stau mass matrix is 
too large and causes one of the stau mass eigenstates to be run below the sneutrino mass. 

Fig. [7| shows the {ttiq, 171112) plane in the z/CMSSM for the two types of models considered 
(corresponding to the normal and inverted hierarchies). As in Figs. [3] and O the brown shaded 
region is excluded because the LSP is the charged partner of the stau. In the left portion of each 
plane, there are areas where one or more of the sparticles are tachyonic. This relatively large region 
is shaded pink. There is also a region which is excluded by measurements of 6 ^ S7 and these 
(shaded green) also exclude low values of mi/2- In the background of each panel, there is a light 
pink shaded region where the value oi — 2 is in agreement with the observed discrepancy to 
within 2a. In the right panel, this region is more prominent and the area within the dashed contour 
satisfies the — 2 constraint within 1 a. The dark blue shaded region in each figure is a direct 
result of the z/CMSSM and corresponds to the area where a sneutrino is the LSP. This region is also 
excluded [47]. Finally, there is the region where the relic density matches the WMAP determination 
as shown by the thin turquoise strips which track either the stau or sneutrino LSP regions. The 
relic density in each case is brought to an acceptable value through neutralino coannihilations with 
staus [48] and/or sneutrinos [49]. 

For our normal hierarchy model shown in the left panel of Fig. [TJ the stau coannihilation region 
is adjacent to the tau sneutrino NLSP regions (the thin blue line shows the position of the contour 
mc,/mf = 1). Analogous to the r sneutrino coannihilation regions, e//i sneutrino coannihilation 
regions appear in our second toy model shown in the right panel of Fig. [71 In this case, the effects 
of the right-handed neutrinos in reducing the left-handed slepton mass are more prominent than in 
the model with a normal neutrino mass hierarchy. For instance, sneutrino coannihilation regions 
show up for Aq > TOOGeV for our inverted mass hierarchy toy model with M^-^ = M^^ = lO^^GeV, 
rriiy^ = TTiy^ = O.OSeV, while Aq > lOOOGeV is required for our normal mass hierarchy toy model 
with M7V3 = lO^^GeV, m^^ = O.OSeV (and other parameter values are fixed as in Fig. [7]). For 
the electron/muon sneutrino coannihilation regions, the WMAP strip transition corresponds to 
the abrupt change from the mostly horizontal right-handed stau NLSP region to the more vertical 
left-handed electron/muon sneutrino NLSP region. For the tau sneutrino NLSP, in contrast, the 
transition from the stau coannihilation region to the tau sneutrino coannihilation region corresponds 
to the smooth change of the stau chirality which is initially right-chiral dominated with an increasing 
left-chiral component for smaller mi/2, before it eventually shifts to the tau sneutrino NLSP region. 
This is expected as the left-handed slepton doublet becomes light due to a large neutrino Yukawa 
coupling in a region with small mi/2- The tau sneutrino typically accompanies an almost degenerate 
left-handed stau state with their small mass splitting coming from D-term contributions, hence 
resulting in nearly parallel contours for rriy^/m^ = 1 and mf/m^ = 1 in the left panel in Fig. [TJ 
while a more abrupt transition from the (right-handed) stau coannihilation to the (left-handed) 
electron/muon sneutrino coannihilation regions arises in the right panel. 
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Figure 7: The (mo,mi/2) plane in the z/CMSSM. In both panels, tan/5 = 10 and Aq = 1000 GeV. 
In the left panel, M^^ = 10^^ GeV and corresponds to the case of a normal neutrino mass hierarchy 
while in the right panel, Mat^ = = 10^^ GeV corresponding to an inverted hierarchy. Contours 
and shading are described in the text. 

6 Conclusion 

The CMSSM has been studied extensively as a template to see the general features of a simple 
super symmetric model. Given the necessity for neutrino masses, the CMSSM implicitly assumes 
that the additional operators relevant to neutrino masses are sufficiently small that their presence in 
the running of the supersymmetric particle spectrum can be ignored. In effect, from the viewpoint 
of the seesaw scheme, such an extension of the CMSSM assumes a relatively low value for the 
heavy right-handed neutrino masses. The z/CMSSM is an extension of the CMSSM and explicitly 
includes three right-handed neutrinos in a general type-I seesaw setup while keeping the universal 
boundary conditions as in CMSSM. Here, we have studied how the dark matter abundances are 
affected by presence of right-handed neutrinos whose masses are close to the GUT scale. Though 
we ignored flavor mixing and the CP violating phase in the neutrino mass matrix, our two simple 
models illustrate the intriguing effects of GUT scale right-handed neutrinos. These include the 
change in the focus point scale, shifts in the funnel region and the realization of light left-handed 
slepton doublets. Clearly the GUT-scale seesaw mechanism and its effects on low-energy observables 
deserve further study. 

This work was supported by Michigan Center for Theoretical Physics (KK), DOE grant DE- 
FG02-94ER-40823 and the William I. Fine Theoretical Physics Institute (KAO). 
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